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The frequency dependence of dynamical conductivity of the quasi-one-dimensional structures with hydrogen 
bonds is studied on the basis of pseudospin-electron model. It takes into account the proton-electron inter- 
action, external longitudinal field h, the tunneling hopping of protons, electron transfer and direct interaction 
between protons. The dependences of the electron concentration and mean number of protons at the site on 
temperature and external field are obtained. The phase transition lines from uniform phase into charge ordered 
phase are determined. The dependence of dynamical conductivity on temperature and field h and its changes 
at the phase transitions are obtained. 
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1. Introduction 

The properties of molecular and crystalline structures with hydrogen bonds are mainly determined 
by the character of proton redistribution on the bonds. We investigate the microscopic mechanisms of 
charge transfer in such systems on the basis of the proposed pseudospin-electron model |1, 2] that takes 
into account the correlation between the proton displacement and reconstruction of electron states as 
well as the change of their occupancy. This interaction manifests itself as a cooperative proton-electron 
transfer (PET) in a number of experimental works 1 3-10] and it follows also from the results of quantum- 
chemical calculations fU ElUlill . Quantum chemical methods allow us to examine these charge redistribu- 
tions more in detail. The structural and optical studies of the proton transfer in AT-salicylideneanihne 0| 
[loll show that photochromism and thermochromism in these object arise from a proton transfer that is 
accompanied by a configurational change of electron structure. It was shown that the behaviour of pro- 
ton dynamics is quite consistent with the temperature dependence of visible absorption spectra of this 
crystal. If we could construct a molecular conductor based on this type of molecules, the charge trans- 
port might strongly be modulated by the proton motion. Photoinduced proton-coupled electron transfer 
(PCET) is investigated in a number of works Il4ll7ll as one of the mechanisms of energy transformation 
in biological and chemical systems. The effect of a such proton-electron coupling plays an important role 
in passing a proton through the biological membrane in photosynthesis. The design of an electron-proton 
hybrid system using the elements of one-dimensional metal chains, acceptor (or donor) molecules, and 
interchain H-bonds are proposed |4]. A new molecular function is expected to be produced in this system, 
if the motion of proton is closely correlated with the dynamics of the ID electronic states. A similar effect 
is observed in the halogen (X)-bridge mixed-valence transition-metal (M) complexes (M-X-complexes) (3JJ. 
The M-X-complexes [MA2X]Y2 (M = Pt, Pd or Ni) have a one-dimensional (ID) chain structure and adja- 
cent chains are connected by hydrogen bonds. Here X stands for a bridging halogen ion (X = CI, Br or 
J), A for a ligand molecule (e.g. ethylenediamine, cyclohexenediamine), and Y for a counter anion (e.g. 
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Y = Br~, ClOp. The location of the protons on N-H-Y induce additional electron charges on the ions M and 
at some conditions they form a charge-ordered state (CDW) |3]. It is pointed out that the electron-proton 
coupling is capable of controlling the CDW state 1181 ■ 

Pseudospin-electron model was originally proposed to describe the correlated proton-electron charge 
transfer in a single complex with hydrogen bond |1]. This model was later extended to the description of 
the charge transfer in the above mentioned quasi-one-dimensional structures with hydrogen bonds. We 
examined uniform phases @1- In this work we investigate the phase transition from uniform phase into 
charge-ordered phase in such systems. We study thermodynamic properties and the frequency depen- 
dence of dynamical conductivity and its changes at the phase transitions. 

2. Hamiltonian 

The Hamiltonian of quasi-one-dimensional structures which contains chains with hydrogen bonds 
are written down in the form 0|: 

#eff = EL{( £ -M)«/a(/) + g["/a(/)-«/ + l,a(/)]S^/)} 
/ i,cr 

+ E E ^aW) \<4 a (l)aj ta (l') + tf a (l')a ilT tf)] 

IJ'i.j.a 1 1 

+EE^(/) 

l i 
l i 

-lLLJm,w')Sp)sp'). (i) 

Here, the summation along the chains (indices i, j) and the summation over the chains (indices /, I') is 
performed. Pseudospin operator S; describes the proton position in double potential well on the hydro- 
gen bond. We suppose that the transfer along hydrogen bond is dominant: t - f;(j),i+i(j); ni a is operator 
of electron concentration at i lattice site, a is electron spin, /i is chemical potential of electrons. 

The Hamiltonian includes proton-electron interaction (parameter g), electron transfer (parameter f), 
energy of proton tunneling (parameter Q), asymmetry of the local anharmonic potential (parameter h). 
The last term describes proton-proton interaction. 

Pseudospin-electron interaction leads to the effective interaction between pseudospins (between pro- 
tons in our case) and as it is shown in (liljilll it can cause the appearance of a modulated phase with 
doubling of the initial lattice period and can lead to the corresponding charge modulation. The study of 
this phenomenon is the aim of this paper. In a case of double modulation of the lattice period, the crystal 
is divided into two sublattices. We introduce the following notations: r/ a = <S? }, n a = ni a „), (a = 1,2 

is the sublattice index). In the mean field approximation (MF) and after passing to ^-representation the 
Hamiltonian Q) has a form: 

H MF = Ha + H sp + U, (2) 
Ha= E [ £ -^ + s[Va-Vi))]n k>a>a + £ t kiaia a^ aa a kAa , 

k,a,a k,a,a 

H sp = EE {«SJ«(0 -[h + jtjp-g [n a - np)] SlJl)} 

l i,a 
IN 

u = -Njvm - —gim - n 2 )(77i -772), T=LLhi),nn ■ 

1 1 V ] 
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The electronic part of Hamiltonian l[2) is diagonalized by unitary transformation 

a-k,\,a - &k,i,o cos <p + afc, 2 ,o- simp, 

«fc,2,a = -ak,l,ir sin( P + ^k,2,(TCOS<p, (3) 

cos2t/? = — . sin2<p = — . 



Thus, we obtain: 



#el = E ( £ fc,a-^)«fc,a,a, (4) 



E Ka ^ E+{-\) a ^ gHm-m) 2 + t 2 k . 

The spin part of Hamiltonian is diagonalized by unitary transformation: 

S*(Z) = Sf (Z) cos i// + Sf(Z) sin 
S*(Z) = -S*(Z) sin i/z + S^ZJcosi/', 
cosi^ a = [h + Jfip-gln a -np)]/7L a ,smi[f a =n/A a . 



In this case: 



H S p = -ZL^i a u), (s) 



3. Thermodynamic properties 

Using formulae (2}-(5), we can write the equations for electron concentration n a and the average 
mean of pseudospins rj a in sublattices: 

1 ^ fl + COS2<0 , ,„, ,„-l 

- E -^{l + exp[/3(£ fc , a - M )]} 



2 k.a.a v 

1 - cos2cp 



{l + exp^^-ju)]}" 1 , (6) 



2 

= ^/irgK-^^fM, (7) 

' 2X a UfcrJ 

From all the possible solutions of equations l[6)-l[7) we choose the ones that correspond to the mini- 
mum of grand canonical potential <1> in regime of \i — const or minimum of free energy F = <J> + fiN in 
regime n = const. In the MF-approximation: 

cD = -2fcrEln({l + exp[-/J(B w -/i)]}{l + exp[-/J(B w -/i)]}) 



k 

-kTNln 

2 



4cosh( — — ] coshf — — 1 
\2kT) \2kT) 



+ -NJri 1 r]2--Ng(n 1 -n2){7]i-T]2). (8) 

From the relations l[6)-([7) we obtain the equations for 8n-n\-ri2 and 8r\-r\\-r\2, which can play a 
role of the order parameter for a modulated phase. Using these equations we obtain the condition of the 
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appearance of nonzero solutions for 8n and 8i), and the equation for temperature of the second order 
phase transition to the modulated phase. 



2 



— (a z ) + — {h + ]r]Y\--{a z ) 2 
A kT 11 U 



?/--^L^({l + «P[i9(e-l<lfcl-M)]r 1 -{l + e^[/'(e+l<lfcl-M)]r 1 ) 

2 N k *k V 7 



1 = 0. 



(9) 



Here: 



m + n 2 m+m . z . i 



n= ^_^ t n=--^, (a^tanh^j, A= sjih + jtf + tf. 



At certain conditions, the order of phase transition can change from the second to the first one. The 
phase transition lines of the first and second order from the uniform phase to the phase with double 
modulation are shown in figure [T] (a) for different values of the chemical potential /i, (b) for different 
values of parameter / and for = and O = 0.05 eV cases. Transition point of the first kind defined by 
numerical calculation as the point at which the requirement of minimum for thermodynamic potential 
with changing the parameters of the model is transition from a homogeneous solution n\ = n.2, 771 = 772 to 
the modulated with different from zero Sn and 5rj. The lines of phase transitions (PT) of the second order 
are shown bold and the lines of the first order PT are thin. The splitting of the electron band at phase 
transition is shown in figure[2] The temperature dependences of mean numbers of electrons of sublattice 
n\, ri2, and uniform phase no also, along the phase transition line are shown in figure[3] The temperature 
dependences of mean values of pseudospins 771, 772, 770 are illustrated in figure[3]as well. The temperature 
dependences of 8n{T) and 6rj{T) are illustrated in figure|4] These results are obtained for the following 
values of parameters: g = 0.08 eV, t = 0.05 eV, / = 0, \i = 0, Q = and O = 0.05 eV. The results for / * 0, 
fi + are presented in figure[T] Such a choice of the parameter values corresponds to the ones given in fl]. 
All energy characteristics (T,h,E,g,t,J,n, Q) are in eV units. 




(a) (b) 

Figure 1 . The phase transition lines of the first and second order from the uniform phase to the phase with 
double modulation: (a) for different values of the chemical potential \i: 1, 2, 3, 4 — p. - 0, 0.05, 0.08, 0. 12 eV 
(n = 0, J = 0), (b) for different values of parameter /: 1,2, 3,4 — 7 = 0,0.05,0.1,0.2 eV (n = 0.05 eV, /i = 0). 
Parameters T, h are in eV units. 
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Figure 3. The temperature dependence of mean numbers of electrons of sublattice n\, 112, and uniform 
phase no and temperature dependence of mean values of pseudospins rji, 772, 770 along the phase transi- 
tion line, n = 0, (/j = 0, / = 0). 




Figure 4. The temperature dependence of the values Sn and Sri along the phase transition line, 1 — = 0; 
2 — Q. = 0.05 eV, (yu = 0, / = 0). 
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4. Dynamic conductivity of quasi-one-dimensional structures 
with hydrogen bonds 

Calculation of the dynamic conductivity of the structure which possesses the chains with hydrogen 
bonds was carried out according to Kubo formula I22I1 



a((u, T^jj-J dfexp[i(w + i£)r] J dA(/(f-iM)/(0)), 




(10) 



where j is the current density operator 



d is dipole momentum operator 



/(0) = -Lff,d], 
n 



(11) 



I i I i 

that includes electronic and pseudospin (ionic) part. Here 8 is the distance between equilibrium positions 
of a proton on the bond, 8 « 0.40 A. According to quantum-chemical calculations, the effective charge of 



hydrogen z^ ff is equal to a 0.25e 



eff , 



j(t) = eTi Ht j{0)e-T I Ht . 



(12) 



In the molecular field approximation, the operator of current density is split into electronic and pro- 
ton (pseudo-spin) parts 

j = je + jsp- (13) 
The following expressions are obtained for these composites: 

dE a m 



Je = 



2e v 



k,a 



dk z 



■{-\) a 2E a {k)F{k) 



a k,a ak ' a 



2e 



+ -r-£-F(*0 [E\{k) -E 2 {k)] {a+ 2 a k> i + a^at^), 
n k 



(14) 
(IS) 



Calculation of correlation functions in the expression (10) with the use of the Wick's theorem yields 
the following expressions for a real part of conductivity: 



a — <7p + £j 



sp, 



(16) 



where the electronic part has a form: 



a e {w) = 



4;re z 
Nah 2 



PL 



k,a 



dE a (k) 

dk z 



■{-l) a 2E a {k)F{k) 



7 8{o)) 



Ane 2 
Nah 2 



^F(fc) 2 [£ 2 (fc)-£i(fc)l 



1 



h 



8 \oj + -[E 2 {k)-E l {k)] \+8 \w- -[E 2 (k) - E^k)] 



{i + e PlE a {k)-fi]y 

e P[E 2 (k)-n] _ e P[E! (*:)-//] 
{ 1 + e P [ft W) -iA } { 1 + e /3 [Eiik) -m] } 
1 



h 



(17) 



Here 



2 g 2 {m-mf + tl 
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Figure 5. Frequency dependence of the electronic 
part of conductance: 1 — before phase transition 
(uniform phase), 2, 3 — after phase transition 
(modulated structure) along the phase transition 
line; 2 — Q. = 0.05 eV, 3 — Q. = 0; \i = 0, / = 0. 



Figure 6. Frequency dependence of the elec- 
tronic part of conductance with different values 
of proton-proton interaction (parameter /), h = 0, 
H = o, T = 0.03016 (350 K): / = 0;0.05;0.14; 1 — 
n = 0; 2 — £2 = 0.05 eV. 



For the protonic part of conductivity we obtain 

\2 1 l-e"' 31 " 



71 I S re Y 1 l-e~P a 

2a\2h n )a^al+e~P Aa 
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Frequency dependence of the electronic part 
tion line is shown in figure [5] curve 1 — before 
phase transition (modulated structure). At 
the phase transition from uniform to modu- 
lated structure, the conductivity cr e (0), when 
co - 0, is abruptly reduced by two to three or- 
ders of magnitude at low temperatures and 
with increasing temperature the value of the 
jump decreases. Electronic conductivity has 
one peak (at co = 0) in uniform phase, (one 
electronic band is present). We observed the 
splitting of the electron band in a modu- 
lated phase, and electronic conductivity has 
a broad maximum in the frequency region 
co = ^[Ezik) -E\[k)\ as well as a peak in 
co — 0. This broad maximum is placed in 
the lower frequency regions for structures 
with the high proton tunneling frequency 
and stronger direct interaction between pro- 
tons. The static conductivity cr e (0) in a mod- 
ulated phase increases with temperature. 
The dynamical conductivity a e {co) decreases 
with increasing temperature and its maxi- 
mum shifts to a lower frequency region and 
vanishes at the critical temperature when 
there is no modulation. In a homogeneous 
phase, only one peak remains at co — 0. The 
change of the frequency dependence of the 
dynamical conductivity with the parameters 
/ and li is shown in figure [6] and figure [7] 
(here we consider the case h — ). Maximum 
of the dynamical conductivity decreases and 
shifts to lower frequency region with an in- 
crease of / and li. Critical values of these pa- 
rameters exist (see figure [1) when the mod- 
ulated phase vanishes and there remains 
only a peak at co = 0. The temperature de- 
pendence of electronic conductivity cr e (0) is 
shown in figure[8] This part of the conductiv- 
ity is higher for systems with larger proton 
tunneling frequency and it increases with an 
increase of the parameters /. Conductivity 
value is presented in relative units. 

Proton dynamic conductivity has peaks 
at frequencies <w ; - = corresponding to pro- 
tons energies A; on hydrogen bonds. One 
peak (Ao) exists in case of homogeneous 
phase and two peaks (Ai, A2) are present for 
the case of a modulated structure. The de- 
pendence of the energy A, on temperature 
and asymmetry field h along the phase tran- 
sition line is shown in figure [9] 



of the dynamical conductivity along the phase transi- 
phase transition (uniform phase), curve 2, 3 — after 
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Figure 7. Frequency dependence of the electronic part of 
conductance with different values of chemical potential 
fi;h = 0,J=0,T = 0.03016 (350 K): Ll = 0.05; 0.07; 0.078; 
1 — n = 0;2 — n = 0.05 eV. 
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Figure 8. Temperature dependence of the electronic part of conductance of the modulated structure 
<7e(0); n = 0;h = 0;l,2 — / = 0;3 — / = 0.05; 4 — / = 0.12; 1 — O = 0; 2, 3, 4 — O = 0.05 eV. 
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Figure 9. The dependence of the peak-frequencies of the proton dynamical conductivity on temperature 
and longitudinal field along the phase transition line; A ; - = hwf, n = 0; / = 0; (a) — O = 0; (b) — O = 0.05 eV. 
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5. Conclusions 

The possibility of the first or the second order transitions from uniform phase into phase with doubled 
lattice period in the quasi-one-dimensional structures with hydrogen bonds is studied in the framework 
of the proposed pseudospin-electron model. It was shown that pseudospin-electron (proton-electron) in- 
teraction may cause the appearance of charge ordered phase in the structures with hydrogen bonds. The 
electron spectrum is calculated. The dependences of the splitting of the electron spectrum on temperature 
and asymmetry field are investigated. The dependences of the electron concentration and mean number 
of protons at the site on temperature and asymmetry field were obtained. It was shown that abrupt 
changes of these characteristics at the first-order transitions are smaller for the structures with high pro- 
ton tunneling frequency and stronger direct interaction between protons. The phase transition lines from 
uniform phase into charge ordered phase are determined. The dependences of the dynamical conductiv- 
ity on temperature and external field and its changes at the phase transitions are obtained. At the phase 
transition from uniform to modulated structure the static conductivity cr e (0) is abruptly reduced by two 
to three orders of magnitude at low temperatures and with increasing temperature the value of the jump 
decreases. Electronic conductivity has one peak at w — in a uniform phase. In modulated phase, the dy- 
namical electronic conductivity has a broad maximum as well as a peak at w = 0. This broad maximum 
is placed at lower frequencies for the structures with high proton tunneling frequency and stronger di- 
rect interaction between protons. It was shown that the frequency dependence of the proton dynamical 
conductivity has one peak in a uniform phase and two peaks in the charge modulated phase. The model 
can be applied to a description of quasi-one-dimensional structures, the so-called halogen-bridge mixed- 
valence transition-metal complexes |3] in which there are charge modulated states. 
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Oa30Bi nepexofli/i i flu Ha m i h h i B/iacTi/iBOCTi 

KBa3iOflHOBMMipHMX CTpyKTyp 3 BOflHeBHMH 3B'fl3l<aMI/l 

P.fl. dei4iB 

iHCTHTyT cf)i3HKH KOHfleHCOBaHWX CWCTeM HAH YKpaiHH, By/1. CBEHL(il4bKOrO, 1, 7901 1 /lbBiB 

Ha ocHOBi nceBflocniH-e^eKTpoHHOi MO,qe/ii floc/iiflxeHO nacTOTHy 3a/iexHicTb flUHaMHHoT npoBiflHOdi KBa3io- 

flHOBWMipHHX CUCTeM 3 BOflHeBUMW 3B'fl3KaMH. B MOfle/li BpaXOBaHO npOTOH-e/ieKTpOHHy B3aeMOfliK), 30BHiLUH£ 

no3flOBXHe no/ie h, TyHe/iK>BaHHfl npoTOHiB, e/ieKTpoHHe nepeHeceHHfl i npaMy npoTOH-npoTOHHy B3aeMOfliK). 
OipniviaHO 3a^exHicTb e/ieKTpoHHOi KOHLieHTpai4iy i cepeflHboT 3ace/ieHor_Ti npoTOHHux no3ni^iw Bifl TeMnepa- 
Typn i no/in h. OipuMaHO /liHiro cf>a30Bnx nepexofliB 3 oflHopiflHOi c(>a3w pp (f>a3n 3 MOfly/inijieK) 3ap^fly. flooii- 
flxeHO 3a/iexHiCTb flUHaMHHOi npoBiflHOcri Bifl no/in h i TeMnepaTypw Ta iT 3MiHH npw cf>a30Bnx nepexoflax. 

K/ihOHOBi c/ioBa: nceBflOcniH-eneKTpoHHa MOfle/ib, npoTOH-eneKTpoHHa B3aeMO^in, BOflHeBm 3b'r30k, 

npOBiflHlCTb 
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